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Abstract. We present a detailed analysis of a high 
resolution spectrum of the damped Lya system at 
z abs = 2.8112 toward PKS 0528-250. The absorption red- 
shift is slightly larger than the emission rcdshift of the 
quasar. We estimate the column density of H 2 molecules 
N(H.2) ~ 6xl0 16 cm -2 and the fractional abundance of 
Ha, / = 5.4x 10~ 5 . The excitation temperature derived for 
different transitions suggests that the kinetic temperature 
of the cloud is ~200 K and the density n ~ 1000 cm -3 . 
The cloud therefore has a dimension of ~1 pc along the 
line of sight. Since it obscures the broad- line emission re- 
gion, its transverse dimension should be larger than 10 pc. 

We obtain upper limits on the column densi- 
ties of C I (< 10 12 - 7 cm" 2 ) and CO (< 10 13 2 cm" 2 ; 
iV(CO)/AT(H i) < 7xl0~ 9 ). We suggest that the ratio 
iV(H 2 )/iV(C i) is a useful indicator of the physical con- 
ditions in the absorber. Simple photo-ionization mod- 
els assuming solar relative abundances show that radia- 
tion fields with spectra similar to typical AGNs or star- 
bursts are unable to reproduce all the constraints and 
in particular the surprisingly small 7V(C i)/iV(H2) and 
7V(Mg i)/7V(H2) ratios. In view of the models we explored, 
the most likely ionizing spectrum is a composite of a UV- 
"big bump" possibly produced by a local starburst and a 
power-law spectrum from the QSO that provides the X- 
rays. Dust is needed to explain the production of molecules 
in the cloud. The amount of dust is broadly consistent with 
the [Cr/Zn] abundance determination. 
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1. Introduction 

QSO absorption line systems probe the baryonic mat- 
ter over most of the history of the Universe (0 < z Ss 5). 
The so-called damped Lya (hereafter DLA) systems 
are characterized by a very large H I column density 
(iV(H i) > 2xl0 20 cm -2 ), similar to the one usually seen 
through local spiral disks. The case for these systems to be 
produced by proto-galactic disks is supported by the fact 
that the cosmological density of gas associated with these 
systems is of the same order of magnitude as the cosmolog- 
ical density of stars at present epochs (Wolfe 1996). More- 
over the presence of heavy elements (Z ~ 0.1 Zq) and the 
redshift evolution of mctallicity suggest the ongoing star 
formation activities in these systems (Pettini et al. 1997), 
while strong metal line systems have been demonstrated 
to be associated with galaxies at low and intermediate z 
(e.g. Bergeron & Boisse 1991). It has also been shown that 
the profiles of the lines arising in the neutral gas show evi- 
dence for rotation (Wolfe 1996, Prochaska k Wolfe 1997). 
Whether these arguments are enough to demonstrate that 
DLA systems arise in large disks is a matter of debate 
however. Indeed simulations have shown that the progen- 
itors of present day disks of galaxies could look like an 
aggregate of well separated dense clumps at high redshift. 
The kinematics could be explained by relative motions of 
the clumps with very little rotation (Haehnelt et al. 1997, 
Ledoux et al. 1998). Moreover, using HSThigh spatial res- 
olution images of the field of seven quasars whose spectra 
contain DLA lines at intermediate redshifts (0.4 ^ z <^ 1), 
Le Brun et al. (1996) show that, in all cases, at least 
one galaxy candidate is present within 4 arcsec from the 
quasar. There is no dominant morphological type in their 
sample: three candidates are spiral galaxies, three are com- 
pact objects and two are amorphous low surface bright- 
ness galaxies. Therefore, although the nature of the DLA 
systems is unclear they trace the densest regions of the 
Universe where star formation occurs. 
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It is thus surprising that despite intensive searches, 
the amount of H2 molecules seems quite low in DLA 
systems in contrast to what is observed in our own 
galaxy. Two detections of H 2 molecules in high red- 
shift DLA systems have been reported. Recently Ge & 
Bcchtold (1997) have found strong absorptions in the 
^abs = 1-9731 DLA system toward Q 0013-004. They de- 
rive N(R 2 ) = 6.9x 10 19 cm" 2 , 6 = 15 km s-\T ex - 70 K 
and n(H) ~ 300 cm -3 for a total hydrogen column den- 
sity N(H) = 6.4xl0 20 cm -2 . This system has by far the 
largest H 2 abundance / = 2iV(H 2 )/[2iV(II 2 ) + N(U 1)] 
~ 0.22±0.05 observed in high z DLA systems. However 
the exact number should be confirmed using a higher res- 
olution data. Other searches have led to much smaller val- 
ues or upper limits (/ < 10~ 6 , Black et al. 1987, Chaffee 
et al. 1988, Levshakov et al. 1992). Levshakov & Var- 
shalovich (1985) suggested that H 2 molecules could be 
present in the z a b s = 2.8112 system toward PKS 0528- 
250. This claim has been confirmed by Foltz et al. (1988) 
using a 1 A resolution spectrum. The latter authors de- 
rive A(H 2 ) = 10 18 cm" 2 , b = 5 km s"\ T ox = 100 K and 
log N(U 1) = 21.1±0.3. By fitting the damped absorp- 
tion together with the Lya emission from the quasar, 
M0ller & Warren (1993) find log N(H 1) = 21.35. Three 
Lya emission-line objects have been detected within 
100ft-- 1 kpc from the quasar by M0ller & Warren (1993) 
and confirmed by Warren & M0ller (1996) to have red- 
shifts within 200 km s _1 from the redshift of the DLA 
system (z a b s = 2.8112 as measured on the Ni 11 lines by 
Meyer & York 1987). The widths of the Lya emission lines 
are very large (> 600 km s _1 ) and continuum emission 
could be present (Warren & M0ller 1996); this suggests 
that the gas is not predominantly ionized by the quasar 
and that star- formation may occur in the clouds, a con- 
clusion reached as well by Ge et al. (1997). The proximity 
of the quasar makes the case difficult however and careful 
analysis is needed. 

In this paper we use much higher spectral resolution 
data to reanalyze the molecular lines in this system. We 
present the observations in Section 2, the results in Sec- 
tion 3 and discuss the low N(C i)/iV(H 2 ) ratio inferred in 
Section 4. 

2. Observations 

The observations were carried out at the F/8 Cassegrain 
focus of the 3.6 m telescope at La Silla, ESO Chile. 
The spectra were obtained with the ESO echelle spectro- 
graph (CASPEC) and long camera. Three exposures of 
5400 s each were obtained under good seeing conditions 
in November 1995. A 300 line mm -1 cross disperser was 
used in combination with a 31.6 line mm -1 echelle grat- 
ing. The detector was a Tektronix CCD with 1124x1024 
pixels of 15 /im square and a read-out noise of 2.5 elec- 
trons. For each exposure on the object, flat field images 
and wavelength comparison Thorium- Argon spectra were 



recorded. The slit width was 1.6 arcsec corresponding to 
a spectral resolution of R <~ 36000 over the wavelength 
range 3650 - 4850 A. The spectra were binned in the di- 
rection perpendicular to the dispersion axis. The accuracy 
in the wavelength calibration measured on the calibrated 
Thorium-Argon spectra is about 0.03 A. 

The data were reduced using the echelle reduction 
package provided by IRAF. The cosmic-ray events have 
been removed in the regions between object spectra before 
extraction of the object. The exposures were co-added 
to increase the signal to noise ratio. During this merging 
procedure the cosmic-ray events affecting the object pix- 
els were recognized and eliminated. The background sky 
spectrum was difficult to extract separately due to the 
small spacing between the orders in the blue. Instead, we 
have carefully fitted the zero level to the bottom of the nu- 
merous saturated lines in the Lya forest. The uncertainty 
on the determination can be estimated to be 5%. 

We identify all the absorption features with equivalent 
widths larger than 5xFWHMxcr where a is the noise rms 
in the adjacent continuum. We use the line list to identify 
the metal line systems in the spectra and more specifically 
all the absorptions associated with the z a b s = 2.8112 DLA 
system. 

3. Results 

3.1. Low ionization metal lines 



Table 1. Heavy element lines from the damped system 



ION 


log N 


[Z/H] 


Ref 


H 1 


21.35±0.10 
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C 1 


< 12.77 




3 


Ar I 


14.46±0.23 


-1.43 ±0.63 


3 


Mg 1 


<12.8 




2,3 


Mg n 


< 15.88 


< -1.04 


2 


Si 11 


16.00±0.04 


-0.90 ±0.11 


2 


S 11 


15.27±0.06 


-1.03 ±0.04 


2 


Cr II 


13.65±0.12 


-1.38 ±0.16 


2 


Fe 11 


15.46±0.09 


-1.41 ± 0.05 


3 


Ni 11 


13.89±0.03 


-1.71 ± 0.10 


2 


Zn 11 


13.09±0.07 


-0.91 ±0.12 


2 



1 M0ller & Warren 1993; ' z Lu et al. 1996; a This work 



Lu et al. (1996) have obtained column densities of most of 
the elements in the DLA system considered here. In addi- 
tion to the lines identified by Lu et al (1996) we find a pos- 
sible Ar 1 doublet from this molecular cloud (see Fig. 1). 
Since our spectrum extends further into the blue, we could 
also obtain a better estimate of the Fe 11 column density 
using the lines from very weak transitions (e.g. Fe nA1125, 
Fe IIA1127, Fe nA1142 etc., see Fig. 1). We also obtained 
an upper limit to the C 1 column density which will be 
useful for constraining the model parameters. The column 
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z abs - 2 8113 system in Q0538-350 
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Fig. 1. Fit result for few weak heavy element transitions in 
the Zabs = 2.8112 damped system 



densities of different species are given in Table 1. Since in 
the neutral phase of DLA systems the dominant ioniza- 
tion state of the elements listed in Table 1 is either the 
neutral (for argon) or the singly ionized (for the others) 
state, we can estimate element abundances. The values 
we obtain for the S n column density is about a factor of 
two less than that of Lu et al. The reason is that we de- 
tect and fit only the strongest component whereas Lu et 
al. integrate the apparent optical depth from v = -80 to 
300 km s~ x . For consistency of the discussion, we adopt 
the value from Lu et al. Subsequently, since singly ionized 
species can be found well outside the region where hydro- 
gen is neutral, it must be noticed that the values of the 
ratios N(X + )/ N(K i) are upper limits. 

An interesting limit on the Mg I column density can 
be obtained from the Lu et al. data. Mg IA2026 is blended 
with Zn nA2026. However from Fig. 9 of Lu et al., we can 
see that the optical depth of the line is smaller than 0.1. 
If we consider, as for Zn II and Cr II that the absorption 
is spread over 40 km s _1 we can infer that w T < 0.027 A. 
Using the oscillator strength given by Morton(1991), we 
find N(Mg i) < 6.6 xlO 12 cm" 2 . 

Formation of H2 molecules is most efficient in the 
presence of dust grains (see Section 3.2). In DLA sys- 
tems the amount of dust is usually estimated from the 
ratio AT(Zn n)/AT(Cr 11) assuming that, like in the ISM of 
our galaxy, Zn is not heavily depleted onto dust grains 
whereas Cr is depleted (Meyer et al. 1989; Pettini et al. 
1994, 1997; Lu et al. 1996). Some controversy has arisen 



however about the presence of dust in DLA systems. Lu 
et al. (1996) have argued that nucleosynthesis alone can 
explain the element abundance ratios observed in DLA 
systems and that the presence of dust is thus question- 
able whereas Pettini et al. (1997) have claimed that even 
though nucleosynthesis history of the gas must play a role, 
the consistent depletion levels of the refractory elements 
indicates the presence of dust. We discuss the abundance 
ratios of Table 1 considering the two possibilities in turn. 

It is known that Zn and S are not heavily depleted 
into dust grains (Sembach & Savage 1996) and should re- 
flect the gas-phase abundances which is thus of the order 
of 0.1-0.2 Zq. The average value of [Cr/Zn] in the diffuse 
interstellar clouds with molecular fraction log /(H2) be- 
tween -6 and -4 is -1.18±0.53 (Roth & Blades 1995). The 
observed [Cr/Zn] = -0.47 in the z = 2.8112 system sug- 
gests that about 70% of Cr is depleted onto dust grains in 
the damped system compared to 99% in the diffuse inter- 
stellar medium of our galaxy. The observed Ni Il/Cr II in 
the 2.8112 system is consistent with the ratio in the inter- 
stellar medium. Silicon has the same gas-phase abundance 
as zinc (see Table 1) which is consistent with the much 
lower depletion of silicon as compared to iron, chromium 
and nickel in the interstellar medium. The non detection 
of Mg 1 is intriguing. Although the upper limit on Mg 11 
is consistent with the other element abundances, some de- 
pletion could be needed (see Section 4). All these are con- 
sistent with the presence of dust in the system with small 
departure of the total element abundance ratios from the 
solar values. 

It is interesting to note however that in our galaxy, 
stars with [Fe/H] ~ -1, have [Si/Fe] ~ [S/Fe] ~ 0.4 (e.g. 
Frangois 1987) and [Cr/Fe] ~ [Ni/Fe] - (Magain 1989, 
Gratton & Sneden 1991) which is approximately what 
we observe. Thus the assumption of depletion onto dust 
grains may be questionable. The upper limit on mag- 
nesium is consistent with [Mg/Fe] ~ 0.4 (Magain 1989). 
However a better limit on this element would definitively 
be very helpful. However if we accept the gas-phase metal- 
licities as the true metallicities, then the observed value 
[Zn/Fe] = 0.5 is not consistent with the value observed in 
the same stars ([Zn/Fe] = 0, Sneden & Crocker 1988). 

In the case of the z a b s = 2.8112 system toward 
PKS 0528-250, we can investigate further this issue by 
asking the question whether dust is needed to produce 
the observed amount of H 2 molecules (see Section 4). 

3.2. H2 molecules 

The consequence of the spectral resolution being nine 
times larger than the previous published data can be ap- 
preciated by comparing Fig. 2 with Fig. 4 of Foltz et al. 
(1988). In our data, a large number of molecular lines are 
left unblended making the determination of the line pa- 
rameters much reliable. 



We have derived the iV(H 2 ) and b values using a Voigt 
profile fitting code (Khare et al. 1997). The oscillator 
strengths of different transitions given in Morton & Diner- 
stein (1976) have been used. Only the clean and unblended 
lines from each rotational level are fitted. The best fit is 
obtained with a single component model although a few 
discrepancies between the fit and the data may argue for 
a two-component model (see also the comment by Jenkins 
& Peimbert f997 about this system). The resulting fit is 
shown in Fig. 2 for a few transitions. The fit parameters 
are given in Table 2 where the second, third and fourth 
columns give the statistical weights of different levels, the 
column density and the number of lines used in the model 
respectively. The best value for the velocity dispersion is 
b = 9.1 km s _1 . 



Table 2. H2 column densities and excitation temperatures 



LEVEL 


g 


Nj 

(10 16 cm" 2 ) 


Number 


T 

1 ex 

(K) 


J cx 


J = 


1 


0.965±0.330 


4 






J = 1 


9 


4.110±0.009 


7 


229±95 


0-1 


J = 2 


5 


0.417±0.063 


8 


208±36 


0-2 


J = 3 


21 


0.216±0.029 


3 


225±08 
224±09 


0- 3 

1- 3 


J = 4 


9 


< 0.03 




< 300 


0-4 



Table 2, the excitation temperatures derived for different 
levels are similar, suggesting that T cx indeed equals the 
kinetic temperature which is therefore close to ~ 220 K in 
this molecular cloud. 

H2 MOLECULE IN THE DAMPED SYSTEM at z = 3.8112 
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The total column density N(R 2 ) in the cloud is 
6xl0 16 cm -2 and the measured fractional H2 abundance 
is / = 2 AT(H 2 )/[2 7V(H 2 ) + N(U 1)] = 5.4xl0" 5 . This 
value is obtained using N(H 1) = 2.2 xl0 21 cm~ 2 (M0ller 
& Warren 1993), and assuming that all the neutral hy- 
drogen contributing to the damped Ly a absorption orig- 
inates from the same region as H 2 . If this were not the 
case, / would be higher. The molecular fraction obtained 
here is two orders of magnitude less than the values quoted 
by Foltz et al. (1988). However as noted by the latter au- 
thors, the column density is very sensitive to the Dopplcr 
parameter and their value should be considered as an up- 
per limit. 

In the case of thermal equilibrium, the ratio of column 
densities of two rotational levels are given by the Boltzman 
equation, 



N(J) _ gj_ 
N(0) g 6XP 



BJjJ+V) 



(1) 



where gj is the statistical weight of any level J in the 
vibrational ground state, T ex is the excitation temperature 
and the constant B/k is ~ 85 K. The values of T ex derived 
for different transitions are given in column 5 of Table 2. 
Lines from higher "J" levels (i.e. J >4) are weak and we 
could derive only upper limits to column densities. 

It is believed that the J = 1 level is mostly populated 
by collisions so the excitation temperature, T i, is roughly 
equal to the kinetic temperature. As can be seen from 



Fig. 2. Fit results for a few rotational transitions of the H2 
Lyman absorption bands in the 2 abs = 2.8112 damped system 



As suggested by Jura (1975) the column density ratio 
N(2)/N(0) and N(3)/N(l) are controlled by the kinetic 
temperature when n(H)>20 and n(H)>300 cm -3 respec- 
tively Our limit on the population of the J — 4 rotational 
level do not add much constraint on the density. How- 
ever, Gc & Bcchtold (1997) report a new determination 
N(4)/N(0) = 4.5xl0- 3 by Songaila & Cowic (1995). This 
new value implies an excitation temperature T 4 = 224 K. 
The calculations by Nishimura (1968) indicate that all this 
is consistent with a kinetic temperature of T ~ 220 K and 
a density n ~ 1000 cm -3 . 

Using the asymptotic approximation given by 
Jura(1974) for the ratio of the molecular to neutral hy- 
drogen densities, n(H 2 )/n(H), as a function of the optical 
depth within a molecular cloud, we can derive an estimate 
of the photo-dissociation rate of molecules, I, due to the 
incident UV flux, 




N 



2</> 



(2) 



where S = 4.2xl0 5 , n is the hydrogen density, R is the 
rate of molecule formation onto dust grains; N is the total 
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H I column density and <f> is the mean molecular frac- 
tion observed through the cloud. If we scale R with the 
metallicities from the value observed in the ISM, then we 
can estimate that R ~ 3xl0~ 18 cm 3 s _1 . This may well be 
an upper limit since the dust to gas ratio in the damped 
system is smaller than in the ISM (see Secion 3.1). Us- 
ing the density and the molecular faction derived above 
(i.e., n ~ 1000 cm~ 3 and / <~ 5.4xl0~ 5 ) we estimate the 
photodissociation rate of molecules , I <~ 3.2xl0~ 8 s _1 , 
and the photo-absorption rate in the Lyman and Werner 
bands, O ~ 7/0.11 - 3xl0~ 7 s" 1 (Jura 1975). This large 
value results from the large H I column density albeit small 
molecular fraction seen in z a b s = 2.8112 damped system. 
We are confident that /, n and the metallicities are known 
within a factor of two. 

The estimated value of the photo- absorption rate, /3 Q , 
is 3xl0 -7 s _1 . This value is an order of magnitude higher 
than the values obtained by Ge & Bechtold (1997). It 
is also two orders of magnitude higher than the average 
value for the interstellar radiation field in our galaxy (i.e. 
0o = 5xl0~ 10 s _1 ). The average value for the intergalac- 
tic medium radiation field calculated for J 2 i(912 A) = 
1, is P = 2xl0~ 12 s _1 . Thus the radiation field in the 
z = 2.8112 absorbing cloud is not dominated by the in- 
tergalactic UV radiation field. The most probable sources 
are either the radiation from the quasar or from the young 
stars in the molecular cloud itself. 

If we take quasar alone (see however Section 4) to be 
the ionizing source then, using tob = 18.17, q = 0.5, 
H = 75 km s _1 Mpc -1 , the optical to UV spectral index 
a = — 1 (f v oc v~ a ) and the upper limit on /3 Q , we estimate 
that the H2 cloud should be at a distance larger than 
10 kpc from the quasar. 

If we assume that the H 2 lines are thermally broad- 
ened then the velocity dispersion of 9 km s _1 gives a ki- 
netic temperature of ~ 2000 K. This is higher than the 
kinetic temperature inferred from the excitation temper- 
atures. Thus it is possible that there are more than one 
component in the H2 absorption. This may introduce a 
small error in the total H 2 column density. However the 
lines are not heavily saturated and hence the errors cannot 
be large. 

3.3. CO molecules 

In the galactic diffuse clouds, N{CO) ~ 10~ 7 xA(H). As 
the high rcdshift DLA systems have, on average, one tenth 
of the solar metallicity we would expect the N(CO)/N(H) 
ratio to be lower than 10~ 7 . Indeed optical and radio 
searches (e.g. Levshakov et al. 1989; Wilkind & Combes 
1994 and references therein) fail to detect CO both in 
emission and absorption in high z damped absorbers. 

However CO emission is detected from a few high red- 
shift quasars and ultra-luminous infrared galaxies (see e.g. 
Barvainis et al. 1994, Ohta et al. 1996; Omont ct al. 1996; 
Guillotcau et al. 1997). These detections suggest that a 



large amount of mass in host galaxies (or surrounding re- 
gions) of some of the high redshift AGNs is in the form of 
molecules. 

We have searched for the presence of absorption due to 
CO A-X molecular bands in a FWHM — 0.6 A spectrum 
obtained by Petitjean & Bergeron (1994) primarily for 
studying intervening C iv systems (see their paper for 
the observational details). No strong CO absorption is de- 
tected from the z = 2.8112 system (see Fig. 3). 

EXPECTED POSITIONS OF CO A X BANDS IN z = 2-8112 ABSORBER 
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Fig. 3. Plot showing the spectrum in the rest wavelength of 
different A-X bands in CO due the z a b s = 2.8112 damped sys- 
tem. The centres of each band is marked with an arrow. 



In order to improve the signal to noise we followed the 
stacking method prescribed by Levshakov et al. (1989). 
The resulting combined spectrum at the rest frame of the 
2-0 band is given in the last panel of the Fig. 3. No clear 
line is seen at the expected position. Using the weighted 
oscillator strength of all the bands and a velocity disper- 
sion of 9 km s _1 , similar to H 2 , we obtain a 3cr upper 
limit for the column density of AT(CO) ~ 1.5 x 10 13 cm~ 2 . 
This limit is consistent with other determinations (Ge et 
al. 1997). This is a factor of two higher than the upper 
limits obtained for the damped systems along the line of 
sight to PHL 957 and 1331+170 (Levshakov et al. 1989). 
This limit leads to iV(CO)/AT(H)< 7 x 10~ 9 . 

Federman et al. (1980) have estimated the CO column 
density along the lines of sight towards 48 bright stars. 
The results of their analysis suggest that CO is readily 
detected when JV(H 2 ) > 10 19 cm~ 2 . Thus the non detec- 
tion of CO in the z — 2.8112 system is consistent with 
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the -/V(CO)/7V(H 2 ) ratios seen in the diffuse interstellar 
medium. They also noted that N(C i) is, on average, 10 
times larger than N(CO). If N(C i) is greater than N(CO) 
in the system under consideration also, then the upper 
limit on N(C i) suggests that the value of iV(CO) should 
be an order of magnitude less than the derived upper limit. 

4. The low 7V(C i)/iV(H 2 ) ratio 

The ratios N(G i)/AT(H i) < 2.2 1CT 9 and 7V(H 2 )/7V(H i) 
= 2.7 derived from this study are surprisingly low. 
Indeed in the interstellar medium of our Galaxy, all the 
clouds with log N(H i) > 21 have log iV(H 2 ) > 19 and 
log N(C i) > 14 (Jenkins & Shaya 1979). It thus may be 
interesting to investigate the reasons for such low ratios. 

Formation of H2 is expected on the surface of dust 
grains if the gas is cool, dense and mostly neutral, and 
from the formation of negative hydrogen if the gas is warm 
and dust free (see e.g. Jenkins & Peimbert 1997). Given 
the neutral hydrogen column density of this system and 
the evidence for depletion onto dust grains (see however 
Warren & M0ller 1996), the first process should dominate. 
Destruction is mainly due to UV photons. The effective 
photodissociation of H 2 takes place in the energy range 
11.1-13.6 eV, through Lyman-Werner band line absorp- 
tion. Since the ionization potential of C 1 is 11.2 eV, it is 
sensitive to the same photo-ionization field as H2 in the 
neutral phase. If the H 2 lines are saturated then the cloud 
is self-shielded against further photodissocation of H 2 . In 
this case however large values of 7V(H 2 ) to N(C 1) ratio 
are easily produced. 

We run grids of photoionization models using 
CLOUDY (Fcrland 1993). These models take into account 
photo-dissociation of H 2 molecules using the approxima- 
tions outlined by Ticlcns and Hollcnbach (1985). The in- 
cident radiation field is averaged over the wavelength in- 
terval 912-1109 A and is redistributed according to the 
interstellar radiation spectrum given by Habing (1968). 

In our calculations we take the gas phase abundances 
to be 1/10 of solar. The grain composition is assumed 
to be similar to that of the Galactic ISM and we use 
the dust to gas ratio as a free parameter. We perform 
the photoionization calculations for different ionizing ra- 
diation spectra and a range of ionizing parameters, U. 
We consider four possible spectral energy distribution 
for the incident radiation field which are (1) a simple 
power-law of index a, F{v) cx v~ a \ (2) a power-law with 
a "big bump" modeled as a black-body of temperature 
50 000 K, the relative contributions of the power-law spec- 
trum and the black-body is characterized by the index 
a ox = log ,F(2500 A)/F(2 keV); (3) a typical spectrum of 
AGNs (Mathews & Ferland 1987); and (4) a starburst 
(Chariot, private communication). The optical to X-ray 
spectral energy distribution of PKS 0528-250 is not well 
defined. Wilkes ct al. (1994) failed to detect X-ray emis- 
sion from this quasar in their Einstein IPC data. They 
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Fig. 4. H2 and C I column densities versus ionization pa- 
rameter for different ionizing spectra (the SEDs are given in 
Fig. |J). The cloud is modeled as a slab of constant density 
n = 100 cm -3 illuminated on one side by an incident flux, re- 
sulting in an ionization parameter U (ratio of the density of 
ionizing photons to the density of hydrogen). Metallicity is a 
tenth of solar and the amount of dust is taken as a fraction of 
the ISM dust content. 
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Fig. 5. Spectral energy distributions of the incident flux ir- 
radiating the cloud from the models presented in the text. 
The small-dashed curve corresponds to a 1 Gyr old starburst; 
the long-dashed curve is the AGN spectrum described by 
Mathews & Ferland (1987); the solid and dotted lines corre- 
spond to power-law spectra with UV big-bumps, modeled as a 
black-body of temperature 50 000 K. 

derive an upper limit on a ox of —1.34. It is known that 
a typical quasar will have cv ox in the range —1.0 to —2.0 
with a mean of —1.5 (e.g. Wilkes et al. 1994). Thus in our 
model we use either a ox = —1.5 or —2.0. 

The predicted column densities of H 2 and C I as a 
function of the ionizing parameter U for different models 
are given m Fig. | (the SEDs are shown for direct com- 
parison m Fig. |). For all models, log N(R i) = 21.35 and 
nn = 100 cm" . We have checked that the exact value of 
riH has little effect on the results. 

To reproduce the small N(C i)/N(H i) ratio, the spec- 
trum must be steep to favor ionization of C I. This is 
why a power-law spectrum with a — -2, a "big-bump" 
spectrum with a ox = -2 and the starburst spectrum can 
reproduce the ratio with somewhat different ionizing pa- 
rameters. The "big-bump" spectrum is favored because 
the X-ray tail maintains the temperature of the neutral 
gas above 100 K. The other models imply too low a tem- 
perature compared to what we derived from the H2 fea- 
tures. However all the models fail to reproduce at the same 
time the H2 column density. It is apparent from the figure 
that steeper spectra could meet the constraints. However 
again the temperature would be too low. Moreover, all 
the models predict log iV(Mg 1) > 13.5. This is a factor of 
seven larger than the upper limit we derived from the Lu 



et al. (1996) data. In the framework of these models, the 
only solution we can find is to decrease the abundances 
of carbon by a factor of two and that of magnesium by a 
factor of seven. This could reflect true relative abundances 
or depletion onto dust grains. 

Dust is needed to explain the H2 column density as 
is apparent from Fig. [|. The H process is more efficient 
than the dust process at temperature T > 500 K and when 
n e is not negligible, thus in the partially ionized region. 
However the very large H 1 column density of this damped 
Lya system requires most of the gas to be neutral where 
this process is negligible because of the low temperature 
resulting from the steep spectrum needed to explain the 
lack of C 1 and Mg 1. 

5. Conclusion 

By fitting the different H2 transitions that we detect in 
a high resolution spectrum of PKS 0528-250, we derive 
log N(H 2 ) ~ 16.78 and T cx ~ 200 K that is most certainly 
the true kinetic temperature of the gas. For this temper- 
ature, the relative populations of rotational levels to 4 
indicate that the density is of the order of 1000 cm~ 3 . 
Therefore the dimension of the neutral cloud along the 
line of sight is less than 1 pc. It must be noticed that the 
damped absorber must cover the broad line region. Indeed 
there is no residual flux in the bottom of the Lya absorp- 
tion line which completely absorbs the Lya emission from 
the quasar over more than 5000 km s _1 (M0ller & War- 
ren 1993). Since the dimension of the BLR in QSOs can 
be approximated as R ~ 0.3 Ljjg 5 where L i6 is the bolo- 
metric luminosity in units of 10 46 erg s _1 (Collin, private 
communication). The radius of the BLR in PKS 0528-250 
is thus of the order of 10 pc. The transverse dimension of 
the damped cloud must be thus larger than 10 pc and the 
cloud must be quite flat. 

We can derive an upper limit on the transverse dimen- 
sion of the cloud by interpreting the non-dection of red- 
shifted 21 cm absorption (Carilli et al. 1996) as an effect 
of partial covering factor of the continuum radio emis- 
sion by the cloud. Indeed the size of the radio source is 
1 arcsec or 5h^ kpc. If the spin temperature is equal to 
the kinetic temperature we derived in sect. 3.2, this im- 
plies that the covering factor should be less than 0.3 and 
thus the radius of the cloud along the transverse direction 
is less than 1 kpc. 

We determine upper limits for the C 1, Mg 1 and CO 
column densities (log N(C 1) < 12.7, log iV(Mg 1) < 12.8 
and log -ZV(CO) < 13.2). Based on simple photoioniza- 
tion models we conclude that (i) no simple model can 
reproduce at the same time the low N(Ci)/N(Hi), 
7V(Mg i)/N(R 1) ratios and the presence of molecules at 
the level observed; (ii) steep spectra are required to re- 
produce the low N(C i)/N(H 1) and 7V(Mg i)/N(H 1) ra- 
tios but they predict temperatures smaller than 100 K in 
conflict with the excitation temperature derived from the 
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H2 transitions; (iii) the only way to keep the temperature 
as high as 200 K is to allow for some X-ray flux heating 
the gas; (iv) in the framework of these models, dust is 
needed to produce the observed amount of molecules; (v) 
the gas phase abundances of carbon and especially mag- 
nesium should be smaller than 0.1 of solar. 

In view of the models we explored, the most likely 
ionizing spectrum is a composite of a UV-"big bump" 
possibly produced by a local starburst and a power-law 
spectrum from the QSO that provides the X-rays. 
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